Nonlinear random resistor diode networks and fractal dimensions of directed 

percolation clusters 

Olaf Stenull and Hans-Karl Janssen 
Institut fur Theoretische Physik III 
Heinnch-Heine- Universitat 

Universitatsstrafle 1 
40225 Diisseldorf, Germany 
(February 1, 2008) 

We study nonlinear random resistor diode networks at the transition from the non percolating 
to the directed percolating phase. The resistor-like bonds and the diode-like bonds under forward 
bias voltage obey a generalized Ohm's law, V ~ I r . Based on general grounds as symmetries and 
relevance we develop a field theoretic model. We focus on the average two-port resistance, which 
is governed at the transition by the resistance exponent cp r . By employing renormalization group 
methods we calculate <f) r for arbitrary r to one-loop order. Then we address the fractal dimensions 
characterizing directed percolation clusters. Via considering distinct values of the nonlinearity r, we 
determine the dimension of the red bonds, the chemical path and the backbone to two-loop order. 
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I. INTRODUCTION 

Percolation ]jj describes the passage of an influence through a medium which is irregularly structured in the sense 
that the influence can propagate through some regions whereas it cannot pass other areas. Prominent examples for 
such media are computer networks like the internet where information propagates and irregularity can be caused by 
random switch failures or other technical problems. 

A particularly simple percolation model is the random resistor network (RRN). In this model the irregular medium 
is given by a, say hypercubic, lattice in which bonds between nearest neighboring sites are randomly occupied with a 
probability p. The influence corresponds to an external current /, which is injected at a terminal site x and withdrawn 
at another terminal site x' . Depending on the occupation probability p the resistors (bonds) are likely to either be 
isolated or to form clusters. Two sites belong to the same cluster if they are connected by a path of bonds and hence 
current can flow between them. At low p two infinitely separated terminal sites x and x' are not connected by such 
a path and the network behaves as an insulator. For large p, on the other hand, many paths between x and x' may 
exist and the network is a conductor. Therefore, at some probability in between, a threshold p c must exist where 
for the first time current can percolate from x to x' . The threshold probability is called the percolation threshold. 
Since it separates the conducting and the insulating phase, it is also referred to as the critical probability. In RRNs 
the influence can percolate through occupied bonds in all directions. The resulting clusters are typically isotropic in 
space. This kind of percolation is referred to as isotropic percolation (IP). The linear extend of the isotropic clusters 
can be characterized by the correlation length £ ~ \p — p c | , where v is the correlation length exponent of the IP 
universality class. 

Directed percolation (DP) [g] is an anisotropic variant of percolation. The bonds function as diodes so that the 
current can percolate only along a given distinguished direction. The critical properties of isotropic and directed 
percolation are very different. Typical DP clusters are anisotropic and they are characterized by two different correla- 
tion lengths: £|| (parallel to the distinguished direction) and £j_ (perpendicular to it). As one approaches the critical 
probability, the two correlation lengths diverge with the exponents v\\ and is± of the DP universality class. 

The apparent success of DP might be attributed to the fact that it is perhaps the simplest model resulting in 
branching self-affine objects. It has many potential applications, including fluid flow through porous media under 
gravity, hopping conductivity in a strong electric field B], crack propagation Q, and the propagation of surfaces at 
depinning transitions . 

DP has a dynamic interpretation in which the distinguished direction is viewed as time. A DP cluster then represents 
the history of a stochastic process. In this dynamic interpretation the DP universality class is the generic universality 
class for phase transitions from an active to an absorbing inactive state. For example the epidemic spreading of an 
infectious desease without immunization [|| may be described by DP 0. Moreover, DP is related to self-organized 
critical models |8|. 
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In the early 1980's Redner |§-|TT| introduced the random resistor diode network (RDN) which comprises both, 
IP and DP. A RDN is a bond percolation model where nearest-neighbor sites are connected by a resistor, a positive 
diode (conducting only in the distinguished direction), a negative diode (conducting only opposite to the distinguished 
direction), or an insulator with respective probabilities p, p+, p_, and q = l—p—p+—p-. In the three dimensional phase 
diagram (pictured as a tetrahedron spanned by the four probabilities) one finds a nonpercolating and three percolating 
phases. The percolating phases are isotropic, positively directed, or negatively directed. Between the phases there are 
surfaces of continuous transitions. All four phases meet along a multicritical line, where < r := p+ = p- < 1/2 and 
V = Pcif)- O n the entire multicritical line, i.e., independently of r, one finds the scaling properties of usual isotropic 
percolation (r = 0). For the crossover from IP to DP see, e.g., Ref. ]l2| ]. In this paper we focus exclusively on the 
vicinity of the critical surface separating the non-percolating and the positively directed phase. 

An important notion in the theory of RDN is the average resistance Mr(x, x') between two connected terminal sites 
x = (x±,t) and x' = (x' ±1 t') of the network. The critical behavior of this average resistance is well known [jl3 14|. If 
Mr(x, x') is measured for example in the distinguished, timelike direction then it scales like 

M R (x,j)~\t-e\+ ,u * , (i.i) 

with <j) being the DP resistance exponent. Here in this paper we consider a generalized RDN in which the resistor- like 
bonds and the diode- like bonds under forward bias voltage obey a generalized Ohm's law, V ~ I r . Our motivation to 
assume this non linear Ohm's law is twofold. First, real circuit elements have non linear current- voltage characteristics. 
This is obviously the case for diodes but it also applies to a certain extend to resistors, in particular for large currents. 
Our main motivation is, however, that the generalized average resistance Mj^ix, x') is related for certain values of 
r to important fractal substructures of DP clusters. This relation provides us with an elegant way to determine the 
fractal dimensions of the red bonds, the chemical path, and the DP backbone. Parts of this work have been presented 
briefly in Ref. ||. 

The plan of presentation is the following. In Sec. [0] we give background on non linear RDN and set up our field 
theoretic model. Section [nj sketches our renormalization group improved perturbation calculation. We derive the 
scaling behavior of Mr (x,x') which is governed by a generalized resistance exponent <p r . We present our one-loop 
result for <fi r . In Sec. [fv[we calculate the fractal dimensions of the red bonds, the chemical path, and the DP backbone 
by considering the limits r — > oo, r — > + , and r — > — 1 + , respectively. Section [v| contains our conclusions. 

II. THE MODEL 
A. Generalized Ohm's law 

Consider a <i-dimensional hypercubic lattice in which the direction n = l/vo(l, . . . , 1) is distinguished. Assume 
that the bonds ^ between two nearest neighboring sites i and j are directed so that bu j\ ■ n > 0. 



Suppose that the directed bonds obey the non-linear Ohm's law 15 

oib.. , (V b ) V b V b , =I b( . , , (2.1a) 

or, equivalently, 

V b = p b (l b ) h \h T^ 1 ■ (2.1b) 

Vb { . = Vj — Vi, where V% denotes the potential at site i, is the voltage drop over the bond between sites i and j. I b 
denotes the current flowing from j to i. In the following we drop the subscript whenever it is save. The bond 
conductances ab are random variables taking on the values a, a6 (V), aO (— V), and with respective probabilities p, 
p+, p-, and q. a is a positive constant and 9 denotes the Heaviside function. The exponents r and s are describing 
the non linearity with r = s . The bond resistances pb are related to the conductances via <7& = p^ s . Note that the 
diodes are idealized: under forward-bias voltage they behave for r = 1 as "oh mic" resistors whereas they are insulating 
under backward-bias voltage. We point out that the round brackets in Eq. ( |2.l| ) contain the argument of the bond 
conductance or resistance, respectively. It is important to realize that these quantities depend on the voltages or 
currents by means of the step function and that sign (Vb) = sign (TV) . Hence we may write a b (H) = Ub_ (lb) and 

Pb (lb) = Pb (Vb) ■ 

Assume that an external current / is injected at x and withdrawn at x' . It is understood that x and x' are connected. 
The power dissipated on the network is by definition 
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P = I(V X -V X ,) 



(2.2) 



Using Ohm's law it may be expressed entirely in terms of the voltages as 



P = R r , + (x,x')- 1 \V x -V xl \ s+1 =J2 (T b(Vb) \V k \ s+ =P({V}) 



(2.3) 



The sum is taken over all bonds on the cluster connecting x and x' and {V} denotes the corresponding set of voltages. 
Rr.+ ix, x') stands for the macroscopic resistance when / is inserted at x and withdrawn at x' . Similarly one may 
define R rj -(x, x') as the macroscopic resistance when / is inserted at x' and withdrawn at x. The two quantities are 
related by R r , + (x, x') — R r ,-(x' , x). From the power one obtains Kirchhoff's first law 



U) 



U) 



-Ii 



as a consequence of the variation principle 

d 
dVi 



= 



The summation in Eq. ([2~4| ) extends over the nearest neighbors of i and Ii is given by Ii = I (<5j )X — Si_ x i). 
Alternatively to Eq. (p. 3]), the power can be expressed in terms of the currents as 



(2.4) 



(2.5) 



P = Rr,+ (X, X') \I\ r+1 =J2Pb_ (h) \h\ r+ = P {{I}) 



(2.6) 



ir+l 



= 



with {/} denoting the set of currents flowing through the individual bonds. It is understood that pt (lb) \h 
whenever <jb (lb) — 0. Kirchhoff's second law, saying that the voltage drops along closed loops vanish, can be stated 
in terms of the variation principle 



d 



dim 



p({j«},j)=o, 



(2.7) 



i.e., there are no independent loop currents 1^' circulating around a complete set of independent closed loops. 



B. Relation to fractal substructures 



Now we take a short detour and show that the average resistance Mn r (x, x') is related for specific values of r to the 
mass, i.e., the average number of bonds, of important substructures of DP clusters. The arguments below are well 
established for RRN and we simply adapt them here to RDN. We start with the backbone. The (forward) backbone 
between two sites x and x' is defined, in principle, as the union of all current carrying bonds whe n / is injected at x 



and withdrawn at x' . Consider r — > — 1 + . One obtains immediately as a consequence of Eq. (2J3), that 



R-i t+ (x,x') = lira /]pb(lb) 
r— 1+ — ' 



E 



Pb 



with only those bonds carrying non zero current contributing to the sum on the right hand side. Hence 



lim + M Rr (x,x') ~ M B 



(2.9) 



where Mb stands for the mass of the backbone. 

Now we turn to r — * oo and r — > + following the lines of Blumenfeld and Aharony |l6| . On the backbone between 
two sites x and x' one may distinguish between two different substructures: blobs formed by multi-connected bonds 
and singly connected bonds which are referred to as red bonds. Both substructures are contributing to the resistance 
of the backbone 
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blob 



R r ,+{x,x') = ^2 Pb{lb) 



r+1 red 

b 



(2.10) 



where the sums are taken over all bonds belonging to blobs and over all red bonds respectively. Since sites on a blob 
are multi-connected by definition, < |/|, and thus 



blob 



lim } p b (l b 



r+1 







(2.11) 



In conclusion, M Rr (x, x') is related to the mass of the red bonds M ro d via 

lim M Rr {x,x') ~ M red . 



(2.12) 



Consider now the first site x at some end of a blob. An entering current / splits into currents Ii^ x flowing to nearest 
neighbors i with 



\ii,*\ = <Ti,*(y ai -v i ) \v x ~v t \ s . 



(2.13) 



In the limit s — > oo the ratios vanish whenever di tX (V x — Vi) \V X — Vi\ a < (Tj >x (V x — Vj) \V X — Vj\ s . Thus, 

current flows only through the resistor with the largest ai jX (V x — Vi) \V X — Vi\ . This argument may be iterated 
through the entire blob. One identifies either a single self avoiding chain through which / flows, with 



(2.14) 



being the power dissipated on the chain, or several of such chains with identical power. The expression in Eq. (2.14) 
is minimal for minimal J2bPb ^- e -> * ne curren t chooses the shortest path through the blob and one is led to 



lim M Rr (x,x') — Mmin , 
r— >0+ 



(2.15) 



where M m i n stands for the mass of the chemical path. 



C. Replica technique 

Our aim now is to calculate the average resistance between two ports x and x' which is precisely defined by 

M Rr {x,x') = (x + (x,x')R r , + (x,x')) c /(x+(x,x')) c . (2.16) 

(■■■)c denotes the average over all configurations of the diluted lattice. x+( x i x ') ^ s an indicator function that takes 
the value one if x and x' are positively connected, i.e., if / can percolate from x to x 1 , and zero otherwise. Note that 
(x+(x, x'))c — (x—{x',x))c is nothing more than the usual DP correlation function. 

Now we follow an idea by Stephen [ p~7j and its generalization to networks of nonlinear resistors by Harris ]l8| and 
exploit correlation functions 



G 

of 



(x,x',\) =(y x (x)Tp_ x (x')) iep (2.17) 
^ x (x)=e^,(iX-V x ) , A^O, (2.18) 



as generating functions of M Rr . In writing Eqs. (2.17) and ( 2.1 8| ) we switched to D-fold replicated voltages, Vi — > 
Vi = (v^\ • • • , V^ D ^ , and imaginary currents, Aj = Hi — > = ^A,- 1 ^, • • • , X\ D ^ ■ The correlation functions are given 
by 
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G (», A a) = lzr° f n fl < ] -p f - 7TT P « f }) 

\ J a— 1 [ 



iuj ■ 

Y 



where 



p ({y}) = t p ({ via) })=f:^-(vt ] 



and Z is the normalization 



Z = 



a=l 



a=l b 



v, 



(a) 



s+1 



(2.19) 



(2.20) 



exp 



s + 



(2.21) 



The additional power term ^ J^i V? which we have introduced in Eqs. (2.19) and ( p. 21 ) is necessary to give the 
voltage integrals a well defined meaning. Without this term the integrands depend only on voltage differences and 
the integrals are divergent. Physically the new term corresponds to grounding each lattice site by a capacitor of unit 
capacity. The original situation can be re trieved by t aking the limit of vanishing frequency w — ► 0. 

Because the integrations in Eqs. ( 2.19 ) and ( 2.21 ) are not Gaussian, we emplo y th e saddle point method. The 
saddle point equation is nothing more than the variation princi ple stated in Eq. (2.5). Thus, the maximum of the 
integrand is determined by the solution of the circuit equations ( [2.4| ). Provided that the condition |/| r+1 ^> a holds, 
we obtain, up to an unimportant multiplicative constant which goes to one in the limit D — > 0, 



G 



(x,x',\\ 



exp 



— R r ,+ (X,X') 



(2.22) 



where 



D 



A r (X) = {- ixi 



a = l 



r+1 



(2.23) 



Now we may expand G about A r ( A ) = 0, 



G 



A r A 



(x,x',\j = (x+(x,x')) c I 1 + — —M Rr {x,x') + 



(2.24) 



This shows us that G is indeed the desired generating function from which the average resistance may be calculated 
via 



M Rr {x,x')={x+{x,x')) c l -rS r-G (x,x',\) 

d(X(Y) /(r + l)J v ' M A )=° 



(2.25) 



Here we would like to comment on the nature of A. We work near the limit when all the components of A are equal 
and continue to large imaginary values. Accordingly we set |18 



= i\ , 



(2.26) 



with real Ao and and impose the condition Yl^—i S,^ = 0. The saddle point approximation may be justified by 
demanding 



lAol » 1. 



(2.27) 
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Substitution of Eq. (2.26) into the definition of A r leads to 

D 



A r A 



+1 r(r + l) lf 



Thus, one can justify the expansion in Eq. (2.24) by invoking the conditions 

and 



|A | r+1 « D- 1 



|A | r - 1 | !l «l 



(2.28) 



(2.29) 



Note that the replica limit D — > allows for a simultaneous fulfillment of the conditions ( 2.27 ) and ( 2.29 ). However, 
we will not only rely exclusively on these conditions on A. We will provide several consistency checks for the validity 
of Harris' saddle point approach as we go along and reproduce known results. 



D. Field theoretic Hamiltonian 



Since infinite voltage drops between different clusters may occur, it is not guaranteed that Z stays finite, i.e., the 
limit lim£)—,o Z D is not well defined. This problem can be regularized by switching to voltage variables $ taking 
discrete values on a D-dimensional torus . The voltages are discretized by setting $ = Adk, where A# = $m /M 
is the gap between successive voltages, z?m is a voltage cutoff, k is a £>-dimensional integer, and M a positive integer. 
The components of k are restricted to —M < few < M and periodic boundary conditions are realized by equating 
_ fc( Q ) m od(2A/). The continuum may be restored by taking &m — > oo and A?9 — > 0. By setting i}m = &oM, 
M = to 2 , and, respectively, A$ = i?o/m, the two limits can be taken simultaneously via to — > oo. Note that the limit 
D — ► has to be taken before any other limit. Since the voltages and A are conjugated variables, A is affected by the 
discretization as well: 



A = AAZ , AAAtf = ir/M , 



(2.30) 



where I is a Z?-dimensional integer taking the same values as k. This choice guarantees that the completeness and 
orthogonality relations 



(2M) 



^A,0mod(2MAA) 



(2.31a) 



and 



(2M) D ^ ' ^) ~ ^.5mod(2MAi?) 



(2.31b) 



do hold. Equation ( [2.31 ) provides us with a Fourier transform between the and A-tori. 

After taki ng care of these regularization issues, we now carry out the average over the diluted lattice configurations 
in Eq. ( [2.19] ). This provides us with the effective Hamiltonian 



if rep = - In ^exp 



s + 



(2.32) 



where 



f D 

^ f i?) = hi I q + p JJ exp 



+ p - n ° xp 

a=l 



S + 



fl+1 



s+1 



■ p+ n ° xp 

a=l 



7TT e ( <,< " > ) 



,' + 1 



(2.33) 



G 



In order to proceed further we recall the choice for A made in Eq. ( [2.26 ). Because A and $ are related by Ohm's law, 
we have to make a consistent choice for 



$ (a) = $o + C (Q) , 

with real #q and <^ Q ) , and where X)a=i — 0- Upon imposing the condition 



|t?o| » 



(a) 



(2.34) 



(2.35) 



for all a we may write 
K (V) 



= In { g + pcxp ( -JT~^ W 



! H?o)+0(tf o )exp 



S+l 



\r)\ 



s+l 



i(# )+6(-ti )exp[ ^-M s+1 

s + l 



(2.36) 



where we have introduced the abbreviation |^| s+1 = X)a=i |$^| S+ ■ After doing a little straightforward algebra and 
by dropping a term 



9 (0 O ) 111 [1 - p - p+] + 6 (-i? ) hi [l-p-p- 
which does not depend on the bond conductances we obtain 



K 



tf) = 9 (d ) K+ (tf) + 9 (-t? ) A-_ (tf 



The K± (tf) in Eq. (|2.38|) are given by 



1 + exp 

1 - p - p ± 



s + l 



(2.37) 



(2.38) 



(2.39) 



Note that these are exponentially decreasing functions in replica space with a decay rate proportional to a . In 

/A • ( 0,, •- i> ) K {,)) . (2.40) 



order to refine H rcp towards a field theoretic Hamiltonian we now expand K (jtj in terms of ip: 



K 



A 1? 



Upon exploiting that 9 (z9 ) = & (Ao) we obtain 



A/0 



i r 

2 



+ 2^( A o)-^(-A )] 



(*)+*-(*) 



X , A - A-_ 



where iT-i- ( A ) stands for the Fourier transform of 



k ± (A) 



(2M) 



(2.41) 



(2.42) 



The Four ier transform can be carried out by switching back to continuous currents and expanding the logarithm in 
Eq. (2.39). The so obtained result has the Taylor expansion 



K± (t\ =t±- uv,± A r (tj + 



(2.43) 



where r± and u> r ,± ~ a 1 are expansion coefficients d epen ding on p and p ± with r±(p,p+,p-) = T^(p,p-,p+) and 
w r.±(p,P+,P-) — Wt,^{PtP-,P+)- Now we insert Eq. ( 2.41 ) into Eq. ( |2.32| ) . We also carry out a gradient expansion 
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in position space. This is justified because only nearest neighbor pairs enter in the power P, i.e., the interaction is 
short ranged not only in replica but also in position space. We find 



-Prep 



EE 

A>0 J i 

~P(Ao) 



K, X 



K_ ( A 



-Ao)] 



A - if _ A 



(») & ■ v + ■ • •] (*) 



(2.44) 



with K± (x \ given by Eq. ( |2.43| ) 



We proceed with the usual coarse graining step and replace the 4>\ (i) by order parameter fields i/j^ (x) which 

inherit the constraint A ^ 0. We model the corresponding field theoretic Hamiltonian H in the spirit of Landau as a 
mesoscopic free energy and introduce the Landau-Ginzburg- Wilson type functional 



r — V - uvA, 



(X) +(0(A o )-0(-Ao))vv] if> x (x) 



A#0 



2 W W ^a+a' W + fv|Vx 



A,A',A+A'^0 



(2.45) 



As usual we have neglected all terms which are irrelevant in the sense of the renormalization group. The parameter 
r is the coarse grained ancestor of r + + t_ . It specifies the "distance" from the critical surface under consideration. 
w r ~ (T _1 is the coarse grained analog of u> r ,+ + The vector v lies in the distinguished direction, v = vn. v 

depends as r and w r on the three probabilities p, p+, and p_. For p + = p_ it vanishes. In the limit iu r -> our 
Hamiltonian TC describes the usual purely geometric DP. Indeed TL leads for w r — > to exactly the same perturbation 
series as obtained in [20| - |22]| . In the limit r — > 1 H describes resistor diode percolation as studied in Refs. Jl|[14|]. In 
the remainder of this paper we drop the regularization term proportional to lu for simplicity. 



III. RENORMALIZATION GROUP ANALYSIS 



Now we are in the position to set up a perturbation calculation. This perturbation calculation can be simplified 
from the onset by manipulating H. in such a way that it takes the form of a dynamic functional [p3|-p5[ . We assume 
that v ^ and introduce new variables by setting 

-1/2 , 



n • x 



vpt , if) = \v\~ 



9=\v\ 1/2 9 



(3-1) 



By substituting Eq. (3.1) into Eq. ( 2.45 ) we obtain 



J= d d 



dt 



A#0 



t - Vi - w r A r A 



>(Ao)-0(-A o ))- 



El V 

6 ^ 

A,A',A+A'#0 



s_ x (x±,t)s_ x , (x_L,t)s x+x , (x_L,i) 



(3.2) 



where d± = d — 1. Note that we have neglected a term containing a second derivative with respect to the "time" t. 
This is justified because this term is less relevant than the one with the first "time" derivative which we kept. 

We proceed with standard methods of field theory |^| . From Eq. (3.2) we gather the diagrammatic elements 
contributing to our perturbation series. The first clement is the vertex —g. Dimensional analysis shows that the 
vertex g is marginal in four transverse dimensions. Hence d = d± + 1 = 5 is the upper critical dimension as it is well 
known for DP. The second diagrammatic element is the Gaussian propagator G t, Xj which is determined by the 
equation of motion 



V 2 



w r A r I A 



[0(Ao)-0(-A o )] 



G 



,t,X 



(3.3) 



For the Fourier transformed G^p,t,Xj of G [xj_,t, A), where p is the momentum conjugate to Xx, one obtains 
readily 
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G (p, t, X) = G+ (p, t, Xj + G_ (p, t, A 
The quantities on the right hand side are given by 

G± (p, t,S)j = 9 (±t) 9 (±A ) exp ^tp (t + p 2 - w r A r (jCj 



For the diagrammatic expansion it is sufficient to keep either G+ ^p,i, A^j or G_ ^p,i, A^. 
G+ (p.t.XY 



(3.4) 
(3.5) 

We choose to keep 



A. Nonlinear resistance of Feynman diagrams 

From the vertex — ~g and the propagator G+ ^p, t, Xj we now assemble the Feynman graphs constituting our dia- 
grammatic expansion. As in our previous work on transport in IP |27] |3l|] these Feynman diagrams have a real-world 
interpretation: they may be viewed as being directed resistor networks themselves. This real-world interpretation has 

basically two roots. The first one is that the principal propagator G+ (p,t, Xj decomposes into two parts: 



G+ p,i,A 



(t) 9 (A ) exp -tp (t + p 2 - w r A r (jCj ) 
-0(t)exp[-ip (r + p 2 )] 6 XS . 



(3.6) 



One of these parts is carrying A's and hence we call it conducting. The other one is not carrying A's and accordingly 
we call it insulating. Equation ( |3 . 6| ) allows for a schematic decomposition of the principal diagrams into sums of 
conducting diagrams consisting of conducting and insulating propagators. In Fig. 1 we list the conducting diagrams 
resulting from the decomposition procedure up to two-loop order. The second root of the real-world interpretation 
is that the replica currents A are conserved in each vertex just as currents are conserved in nodes of real networks. 
Hence we may write for each edge i of a diagram, A^ = Ai ^A, {k}^ , where A is an external current and {/?} denotes a 

complete set of independent loop currents. The A-dependent part of each conducting diagram then takes the form 



exp 



pw r ^2uA r (Xi^j 



(3.7) 



Now it is important to realize that J2i (^i) resembles the structure of a power (cf. Eq. (p^)). Thus, we interpret 
the "time" associated with a conducting propagator as its resistance and write 



(3.7) = exp pw r P r ( A, {k} 



(3.8) 

The real- world interpretation provides for an alternative way of computing the conducting Feynman diagrams. To 
evaluate the sums over independent loop currents, 



E 

{«} 



exp 



pw r P r (X, {k}) 



(3.9) 



we employ the saddle point method under the conditions discussed a t th e end of Sec. II C . Note that the saddle point 



equation is nothing more than the variation principle stated in Eq. (2.7). Thus, solving the saddle point equations is 
equivalent to determining the total resistance R r ({ti}) of a diagram, and the saddle point evaluation of (|3.9|) yields 



exp 



Rr({U})pw r A r (A 



(3.10) 



where we have omitted once more multiplicative factors which go to one for D — * 0. A completion of squares in the 
momenta renders the momentum integrations, which remain to be done to compute the diagrams, straightforward. 
Equally well we can use the saddle point method which is exact here since the momentum dependence is purely 

quadratic. After an expansion for small w r A r (jCj all diagrammatic contributions are of the form 







/(p 2 ,t,A r (X)) =I P (p 2 ,t)+I w (p 2 ,t)pw r A r fy +••■ 

Y[dt z [l + Rr ({<»}) P w r A r (X) +••■] £(p 2 ,i;{ij) 



(3.11) 



D (p 2 , i; {ii}) is a typical integrand as known from the field theory of DP [£0] [22] 



B. Renormalization and scaling 



We proceed with standard techniques of renormalized field theory p6[ . The ultraviolet divergences occurring in the 
diagrams can be absorbed by dimensional regularization. We employ the renormalization scheme 



- s = Z 1 l 2 s , 

T = Z^ZtT , 



P P 

w r — > 



Z-'Zpp 



W r = Z p 1 Z Wr W r 



(j, 



(3.12a) 
(3.12b) 
(3.12c) 



where e = 4 — and /i is the usual inverse length scale. The factor G t = {A-K)~ d± > 2 T(1 + e/2), with T denoting the 
Gamma function, is introduced for convenience. Z, Z T , Z p , and Z u are the usual DP Z factors known to second order 
in e p^,p3,0|. In Rcf. (l^] we determined Z w = Z Wl to second order in e. Here, we calculate Z Wr for arbitrary r to 
order e. This calculation is straightforward because we can determine the total resistance of the one-loop diagrams 
by using simple rules. For example, two nonlinear resistors with resistances t\ and £2 added in series have a total 
resistance R r given by 



Rr(tl,t2) — t\ + t 2 , 



whereas two such resistors in parallel give 



Rr(ti,t2) 8 — t 1 s + < 2 ' 



By exploiting Eq. (3.14) find 



it 

Ye 



O u 2 



(3.13) 
(3.14) 

(3.15) 



Calculating Z Wr for general r to higher loop orders appears to be beyond possibility. The reason is, that conducting 
diagr ams l ike C appe ar. The total resistance of these diagrams cannot be determined by using simple rules like 
Eqs. fl3.13p and ( frj ). Instead, one has to solve the set of nonlinear circuit equations which is hardly feasible in 
closed form. 

Now we set up in a standard fashion the renormalization group equation for our problem. The unrenormalized 
theory has to be independent of the length scale pT 1 introduced by renormalization. In particular, the unrenormalized 
connected TV point correlation functions must be independent of /i, i.e., 



M^-Gat ({ X -L> pt,w r A r (X) j ;f,5 







for all N. Eq. ( 3.16 ) translates via the Wilson functions 



K (u) = /J, 



din- 



d In w r 



dp 



7... (u) 



C P (u) = M 

M^-in Z.. 
op 



dp 

91np 



dp 



(3.16) 

(3.17a) 
(3.17b) 
(3.17c) 



where the bare quantities are kept fix while taking the derivatives, into the Gell-Mann-Low renormalization group 
equation 
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d n d d _ d d N 



x G 



N 



N x x , pt, w r A r (Xj | ; r, u, fi\ 







(3.18) 



The particular form of the Wilson functions can be extracted from the renormalization scheme and the Z factors. 

T he cr itical behavior of the correlation functions is determined by the infrared stable fixed point solutions of 
Eq. (3.18). This fixed point u* is readily extracted from the condition (3(u*) = 0. Then Eq. (3.18) is solved at u* by 
the method of characteristics which gives 



G N (jxi, pt, w r A r (X^j | ; r, u, p 
= P' N/2 G N ({lx±,&pt, l c *w r A r (a) } ; J"V, u*,lp 



(3.19) 



where 7* = 7 (u*), k* — k (u*), C» — C p ( u *) > an d C* — Cr ( u *)- To analyze the scaling behavior of the correlation 
functions completely, the solution ( 3. L9[ ) has to be supplemented by a dimensional analysis: 



G N ({ xx, pt, w r A r (\j X ; r, u, p, 
= p d ^ N ^ 2 G N ^ X± ,p 2 pt,p- 



■A, 



(a) } ; n 2 r, u, 1 



(3.20) 



Equation (3.19) in conjunction with Eq. (3.20) now gives 
Gat ^x±,pt, w r A r rXjj ;T,u,pA 



= l {d ^ )N/2 G N ({/xx, I'pt, r^'^WrAr (a) } ; P^T, «*,/*) 



(3.21) 



77 = 7*, z = 2 + £*, and z/x = are the well known the critical exponents for DP which have been calculated 

previously to second order in e |2l| , |2l||l4|| . These DP expon ents, however, are not sufficient to specify the critical 
behavior of the RDN correlation functions completely. In Eq. (3.21) we introduced the additional nonlinear resistance 
exponent 



Vy_ (2 - C) = 1 + 



3 ■ 2' + 2 



Oi 



(3.22) 



Note tha t <p = <pi is in conformity to order e with our result for the resistance exponent for the usual "ohmic" RDN, 
i.e., Eq. ( |3.22 ) satisfies an important consistency check. 

Since we are primarily interested in the critical behavior of the average two-port resistance, we now take a closer 
look at the two point correlation function G = Gi. Equation (3.21) implies for G at r = that 



G (|XX - XxM - t', W r A r (X)) = l d ±+*G (l\*X ~ Xj. M* (* - , l'^ 1 ^ Wr A r (xj) 



(3.23) 



where we dropped several arguments for notational simplicity. In the following we set x^_ = and t' = 0, once more 
for the sake of simplicity. The choice I — |x^| _1 and a Taylor expansion of the right hand side of Eq. (3.23) lead to 



G \x±\,t,w r A r (A 



Ixxl 1 -'""/! 



t 



IXxl^A,,! 



t 



xx 



X |l + W r A r (S^J 

where the f's are scaling functions. Equally well we can choose / = t~ 1/>z which then leads to 
G(\x x \,t,w r A r (\)) =t(i-"-W Va ^ 



(3.24) 



l + m r A,.(A) t^f u 
with other scaling functions fa and f Wt 2 and where i>u = v±z. 



t 



(3.25) 
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Now we can extr act th e critical beha vior of Mr t . For measurements in the distinguished direction we straightfor- 
wardly exploit Eq. ( |3.25| ) via Eq. ( 2.25 ) and find that 



AI 



Rr 



For all other directions we determine Mr t from Eq. ( 3.24 ). With help of Eq. ( 2.25 ) we find that 

M Rr 



\4> T /U ± n I t 

X -L 



Here it is convenient to choose a common length scale L and to express both, |xj_| a nd t, in terms of it: |xj 
t ~ L z . This choice guarantees that the scaling functions f w \ is constant and Eq. (3.27) simplifies to 



(3.26) 

(3.27) 
L and 



M 



Rr 



(3.28) 



IV. FRACTAL DIMENSIONS OF DP CLUSTERS 

In this section we calculate <j) r for r — » + , r — > oo, and r — > — 1 + to two- loop order. This provides us with the 
fractal dimension of the backbone, the red bonds, and the chemical length respectively. 



A. Red bonds 



For self-affine objects the notion of fractal dimension is less straightforward than for self-similar objects. To 
determine, for example, the fractal dimension of the red bonds in DP one considers a d — 1-dimensional hyper-plane 
with an orientation perpendicular to x«. The cut through the red bonds is a self similar object with the fractal 
dimension 



jcut j -I 

"rod — "red 1 



(4.1) 



where <i rc d is the local fractal dimension [|32| of the red bonds. By virtue of Eqs. (2.12) and (3.27) the mass of the red 
bonds scales as 



\^±\ Mux f w ,i 



\\ 



XL 



(4.2) 



Accordingly the mass of the cut scales like 



M cut 



XL 



&-l/l/_L 



X 1 1 fw , 1 



(4.3) 



By choosing once more |xj_| ~ L and xn ~ L z we find that 



M cut 
red 



f j 4> 00 /v±-z 



(4.4) 



This leads via Eq. (O) to 



1 + 4>oo/v±_ - Z . 



(4.5) 



It remains to compute . To do so w e tak e direct advantage of our view of the Feynman diagrams as being resistor 
networks themselves. As argued in Sec. [IB , blobs do not contribute to the total resistance for r — > oo. In analogy 
only singly connected conducting propagators contribute to the total resistance of a diagram, i.e, 



singly 

Rao ({U}) — ^ U i 



(4.6) 



with the sum being taken exclusively over singly connected conducting propagators. The contribution of a diagram 
to the rcnormalization factor Z,„ takes the form 
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singly 



LW (P 



-*)=/ n^' i) (p^;W)E t 

JO 



(4.7) 



Note that a factor ^ in Eq. ( |4.7| ) corresponds to the insertion (cf. Ref. f26|| ) of is 2 into the ?th edge of the diagram. 
Thus, we can generate Iw (p 2 >^) f° r a given conducting diagram by inserting is 2 into its singly connected conducting 
propagators. This procedure is carried out up to two- loop order, i.e., every conducting propagator in Fig. 1 that does 
not belong to a closed loop gets an insertion. The resulting diagrams are displayed in Fig. 2. 

At this point it is instructive to consider the contributions of the diagrams listed in Fig. 1 to Z T . These can 
be generated by inserting is 2 in conducting as well as in insulating propagators. Again, one obtains the diagrams 
depicted in Fig. 2 with the same fore-factors. Consequently, Z WaB 
The same goes for the corresponding Wilson functions £oo and k. 



and Z T are identical at least up to two-loop order. 
From the definition of 6 r it follows that 



2-C 



= 1 + (e 3 ) 



(4.8) 



Upon inserting Eq. (4.8) into Eq. ( |4.5| ) we obtain the prime result of Sec. IVA, 

drod = 1 + 1/W - Z , 



(4.9) 



holding at least to second order in e. By substituting the results of Ref. [^l|,^2| in Eq. (4J)) we obtain the following 
e-expansion for d rc( i'. 



dred 



73 55 /4 
288 ~ 144 n V 3 



(4.10) 



Note that the scaling relation (4.9) holds rigorously and its validity is not restricted to second order in the e- 
expansion. Coniglio |33] proved for IP that the mass of the red bonds scales as 



M red (0, x) ~ |x| 



l/u 



(4.11) 



For IP this is equivalent to saying that the fractal dimension of the red bonds is d rc d = X/v. Since Coniglio's argum ents 
do not rely on the isotropy of the system, they can be adapted to apply to DP. For the DP problem Eq. ( 4.11 ) has 
to be modified to 



M red (0,x) = \x x \ 1/vx f w ,i 



(4.12) 



This in turn leads again to Eq. (4.9) 



B. Chemical length 



Next we address the fractal dimension of the chemical length. Equation ( |2.15| ) in conjunction with Eq. ( 3.27 ) 
provides us with 

M min = \x ± f o/v± f W)1 (t^y) ■ (4-13) 



By applying the same reasoning as in Sec. [V C we learn that the local fractal dimension of the chemical length is 
given by 

eUn = 1 + 0o/W - * • (4-14) 

In order to calculate cf>o we determine the shortest self-avoiding path of conducting propagators connecting the 
external legs of a diagram. Due to the dynamic structure all of these paths for a given diagram have the same total 
resistance that is nothing more than the total "time" between the external legs. Hence, we can choose any self-avoiding 
path connecting the external legs. We work with the diagrammatic expansion depicted in Fig. 3. 

Minimal subtraction provides us with the renormalization factor 
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,2 r 



Z = 1 



4e 32 e 



7 o 9 , 
e- 3+ 2 ln 3 



(4.15) 



Note that as stated in Eq. (4.15) is iden tical t o the field renormalization Z given to two-loop order in Refs. ]2^j22| , |l4j . 
By virtue of the renormalization scheme ( |3.12 ) we deduce that 



pw = pw 



at least to second order in u. Equation ( Hq ) leads via 

d In (p w ) 



dp 



= o 



to 



C P + Co = o 



From the definitions of <p r and z it follows immediately that 

cj) = u± z + O (e 3 ) 



In conjunction with Eq (2.15) this leads finally to 

dmi„ = l + 0(e 3 ) . 

This result is intuitively plausible because the chemical distance in DP is basically equivalent to the "time" t. 



(4.16) 

(4.17) 

(4.18) 
(4.19) 
(4.20) 



C. Backbone 

We conclude Sec. |v| by studying the backbone dimension. By virtue of Eqs. (2.9) and ( 3.27 ) the mass of the 
backbone scales as 

M B = \x ± \^^f wA (j^j . (4.21) 

Accordingly the local fractal dimension of the backbone is given by 

D B = \ + <t>-i/v±-z . (4.22) 



It remains to compute </>_i. Once more we exploit our real-world interpretation. As argued in Sec. |II Bj , the 
resistance of the backbone between two sites x and x' is given by 

R^, + (x,x')=J2pb(lb) , (4.23) 

b 

with the sum running over all current carrying bonds of the underlying cluster. In analogy, the resistance of a 
conducting Feynman diagram is given by 

cond 

R-i {{U}) = U . ( 4 - 24 ) 

i 

where the sum is extending over all conducting propagators of the diagram. This tells us that the contribution of a 
diagram to Z W l takes the form 

„oo cond 

iw{v\t) = / nd^i>(p a ,t s {*j})y;*< . (4.25) 

Hence we can generate Iw (p 2 i*) 01 anv conducting diagram by inserting |s 2 into its conducting propagators. All 
conducting propagators in Fig. 1 get such an insertion. Then it is a matter of simple counting to see that the individual 
contributions cancel each other. Thus, we find 
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Z w _, = 1 + (e 3 ) 



(4.26) 



As a consequence we obtain 



lim <j) r / u± ~ z — rj 



at least to second order in e. Equation ( 4.27 ) leads by virtue of Eq. ( 4.22 ) to 



D B = 1 - rj = d - — 



(4.27) 



(4.28) 



where = v x{d — 1 + 77) /2 is the DP order parameter exponent known to second order in e |^T|,^2|. From the scaling 
relation Eq. (4.28), which is the main result of this section, the e-expansion of Dg is readily obtained by inserting the 
e-expansion for rj [2l|p2]| : 



D B = l + -l\ 



25 161, [A 

1 In - 

288 144 V 3 



e + O (e 3 



(4.29) 



Equation ( 4.28| ) is in agreement with scaling arguments [|4| yielding that the fractal dimension of the transverse 
cut through a DP cluster with local dimension d/isd/ — 1 = d — 1 — (3/v±. The analogous cut through the backbone 
can be viewed as the intersection of the cut through the cluster and the clusters backward oriented pendant |lO|,[35| . 
Hence, the codimension of the backbone cut is twice the codimension [3/v± of the cluster cut, which leads again to 
Eq. ( p8| ). 

It is interesting to compare the e-expansion result to numerical estimates. We are not aware, however, of any 
simulations in which Db itself was determined. Huber et al. p6| presented numerical results for the scaling exponent 
of the backbone mass when measured in the longitudinal direction. In the following we call this exponent Db- 
Formally one can define Db via Mb ~ t Dt 



From Eqs. (2.£), (3.26), and (4.27) it follows that 



D, 



1 



1 



(4.30) 



at least to second order in e. Crudely evaluating the corresponding e-expansion of Db for small spatial dimensions leads 
to poor quantitative predictions. Therefore it is appropriate to improve the e-expansion by incorporating rigorously 
known features. We carry out a rational approximation which takes into account that obviously Ds(d = 1) = 1. 
Practically this is done by adding an appropriate third order term. By this procedure we obtain the interpolation 
formula 



D B » 1 



1 — -J (0.0833 e + 0.0583 e 2 ) . 



(4.31) 



Evaluation in two dimensions leads to Ds{d = 2) = 1.2 ± 0.15, where the error is based on a subjective estimate. 
This result is, within the errors, in agreement with the numerical result j3(| Ds{d = 2) = 1.30 ± 0.03. Our result 
though appears to be somewhat small. 



V. CONCLUSIONS 



In this paper we studied a nonlinear version of resistor diode percolation where Ohm's law is generalized to V ~ I r ■ 
We investigated the critical behavior of the average two-port resistance Mr t at the transition from the non percolating 
to the directed percolating phase. By employing our real-world interpretation of Feynman diagrams we calculated 
the resistance exponent cj> r for arbitrary r to one-loop order. To our knowledge this is the first time that <p r has been 
determined for the RDN while 4> r is known for RRN, also to one-loop order, since the 1980's Extending either 
of these results to higher loop orders seems to be beyond possibility because not all conducting diagrams appearing 
at higher loop orders can be assembled by simply adding resistors in parallel and in series. For these diagrams one 
has to solve the set of nonlinear Kirchhoff 's equations to obtain their total resistance. In closed form, however, this 
is hardly feasible. 

The relation of Mr^ to the mass of the red bonds, the chemical length, and the backbone, respectively, provided 
us with alternative means to extract the fractal dimensions of these substructures of DP clusters. By computing <j) r 
for r — > 00, r — * + , and r — * — 1 + we determined d re d, d min , and Db to two- loop order. 
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The fractality in DP and IP is qualitatively different. DP clusters arc self affine rather than self similar objects. 
Hence, the notion of fractal dimensions in more subtle for DP than for IP. Moreover, DP has a Markovian character 
which is evident in the dynamic interpretation. This Markovian character provides for scaling relations which do 
not have an analog in IP. The DP backbone dimension for example can be expressed entirely in terms of the usual 
(purely geometric) critical exponents of the DP universality class, Db = d — 2f3/v±. Within the renormalization 
group framework such a scaling relation is typically associated with a Ward identity. The fact that W-\ rcnormalizcs 
trivially to two-loop order is reminiscent of this Ward identity. It is an interesting issue for future work to identify 
the Ward identity and its underlying symmetry. Another consequence of the Markovian character of DP is that the 
fractal dimension of the chemical length is identical to one. This is intuitively plausible since the shortest longitudinal 
path through a DP cluster corresponds to the time in the dynamical interpretation. 

The fractal dimension of the red bonds in DP obeys the scaling relation d rcd = 1 + — z similar to the d lcd = 1/v 
for IP. The Ward identities corresponding to cither of these scaling relations are not known to date. Again, this leaves 
interesting and challenging opportunities for future studies. 



ACKNOWLEDGMENTS 



We acknowledge support by the Sonderforschungsbereich 237 "Unordnung und grofie Fluktuationen" of the Deutsche 
Forschungsgemeinschaft . 



[1] For a review on percolation see, e.g., A. Bunde and S. Havlin, Fractals and Disordered Systems (Springer, Berlin, 1991); 

D. Stauffer and A. Aharony, Introduction to Percolation Theory (Taylor&Francis, London, 1992); B. D. Hughes, Random 

Walks and Random Environments (Clarendon, Oxford, 1995). 
[2] For a review on DP see, e.g., H. Hinrichsen, Adv. in Phys. 49, 815 (2000). 
[3] N. Van Lien and B. I. Shklovskii, Solid State Commun. 38, 99 (1981). 
[4] J. Kertesz and T. Vicsek, J. Phys. C 13, L343 (1980). 

[5] See, e.g., A.-L. Barabasi et al, in Surface Disordering: Growth, Roughening, and Phase Transitions, edited by R. Jullien 
et al. (Nova Science, New York, 1992); S. Havlin et al, in Growth Patterns, edited by J. M. Garcia-Ruiz et al. (Plenum, 
New York, 1993). 

[6] See, e.g., J. D. Murray, Mathematical Biology (Springer, Berlin, 1988). 
[7] P. Grassberger, J. Phys. A 18, L215 (1985). 

[8] K. Sneppen and M. H. Jensen, Phys. Rev. Lett. 70, 3833 (1993); 71, 101 (1993); Phys. Rev. E 49, 919 (1994); P. Bak and 
K. Sneppen, Phys. Rev. Lett. 71, 4083 (1993); M. Paczuski, S. Maslov, and P. Bak, Europhys. Lett. 27, 97 (1994). 

[9] S. Redner, J. Phys. A: Math. Gen. 14, L349 (1981); Phys. Rev. B 25, 3242 (1982). 
[10] S. Redner, in Percolation Structures and Processes, edited by G Deutscher et al. (Adam Hilger, Bristol, 1983). 
[11] RDNs were already contained implicitly in the pioneering work of S. R. Broadbent and J. M. Hammersley on percolation, 

Proc. Philos. Soc. 53, 629 (1957). 
[12] H. K. Janssen and O. Stenull, Phys. Rev. E 62, 3173 (2000). 
[13] H. K. Janssen and O. Stenull, Phys. Rev. E 63, 025103(R) (2001). 
[14] O. Stenull and H. K. Janssen, to appear in J. Stat. Phys. 
[15] See also S. W. Kenkel and J. P. Straley, Phys. Rev. Lett. 49, 767 (1982). 
[16] R. Blumenfeld and A. Aharony, J. Phys. A 18, L443 (1985). 
[17] M. J. Stephen, Phys. Rev. B 17, 4444 (1978). 
[18] A. B. Harris, Phys. Rev. B 35, 5056 (1987). 

[19] See, e.g., A. B. Harris and T. C. Lubensky, Phys. Rev. B 35, 6964 (1987). 
[20] J. L. Cardy and R. L. Sugar, J. Phys. A 13, L423 (1980). 
[21] H. K. Janssen, Z. Phys. B 42, 151 (1981). 
[22] H. K. Janssen, J. Stat. Phys. 103, 801 (2001). 

[23] H. K. Janssen, Z. Phys. B 23, 377 (1976); R. Bausch, H. K. Janssen, and H. Wagner, ibid. 24, 113 (1976); H. K. Janssen, 
in Dynamical Critical Phenomena and Related Topics, edited by C. P. Enz, Lecture Notes in Physics, Vol. 104 (Springer, 
Heidelberg, 1979). 

[24] C. De Dominicis, J. Phys. (Paris) Colloq. 37, C-247 (1976); C. De Dominicis and L. Peliti, Phys. Rev. B 18, 353 (1978). 
[25] H. K. Janssen, in From Phase Transitions to Chaos, edited by G Gyorgyi, I. Kondor, L. Sasvari, and T. Tel (World 
Scientific, Singapore, 1992). 



1G 



[26] See, e.g., D. J. Amit, Field Theory, the Renormalization Group, and Critical Phenomena (World Scientific, Singapore, 

1984); J. Zinn- Justin, Quantum Field Theory and Critical Phenomena (Clarendon, Oxford, 1989). 
[27] O. Stenull, H. K. Janssen, and K. Oerding, Phys. Rev. E 59, 4919 (1999). 
[28] H. K. Janssen, O. Stenull, and K. Oerding, Phys. Rev. E 59, R6239 (1999). 
[29] H. K. Janssen and O. Stenull, Phys. Rev. E 61, 4821 (2000). 

[30] O. Stenull, Renormalized Field Theory of Random Resistor Networks, Ph.D. thesis, Universitat Diisseldorf, (Shaker, 
Aachen, 2000). 

[31] O. Stenull and H. K. Janssen, Europhys. Lett. 51, 539 (2000). 

[32] See, e.g., J. Kertesz, and T. Vicsek, in Fractals in Science, edited by A. Bunde and S. Havlin (Springer, Berlin, 1994). 

[33] A. Coniglio, Phys. Rev. Lett. 46, 250 (1981); J. Phys. A 15, 3829 (1982). 

[34] See, e.g., B. Hede, J. Kertesz, and T. Vicsek, J. Stat. Phys. 64, 829 (1991). 

[35] B. M. Arora, M. Barma, D. Dhar, and M. K. Phani, J. Phys. C 16, 2913 (1983). 

[36] G. Huber, M. H. Jensen, and K. Sneppen, Phys. Rev. E 52, R2133 (1995). 



17 




FIG. 1. Decomposition of the primary two leg diagrams (bold) into conducting diagrams composed of conducting (light) 
and insulating (dashed) propagators to two-loop order. It is important to realize that the conducting diagrams inherit their 
combinatorial factor from their bold diagram. For example, the diagrams A and B introduced below have to be calculated with 
the same combinatorial factor | . 




FIG. 2. Diagrammatic expansion in the limit r — * oo. The listed diagrams including their fore-factors can be obtained from 
the conducting diagrams shown in Fig. 1 in two different ways: first, by inserting |s 2 into all singly connected conducting 
propagators and second, by inserting |s 2 into every conducting and insulating propagator. As a consequence, the renormal- 
ization factors Z W(X1 and Z T are identical. The lines stand for conducting propagators evaluated at zero currents, the solid dots 
for is 2 -insertions. 
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FIG. 3. Diagrammatic expansion in the limit r — > + . The meaning of the symbols is the same as in Fig. 2. 
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